Boundary Friction on Molecular Lubricants: Rolling Mode? 
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A theoretical model is proposed for low temperature friction between two smooth rigid solid 
surfaces separated by lubricant molecules, admitting their deformations and rotations. Appearance 
of different modes of energy dissipation (by "rocking" or "rolling" of lubricants) at slow relative 
displacement of the surfaces is shown to be accompanied by the stick-and-slip features and reveals 
a non-monotonic (mean) friction force vs external load 

PACS numbers: 46.55.+d, 81.40.Pq 



INTRODUCTION 

In the modern tribology a still increasing interest is put 
to the studies of wearless friction on atomically smooth 
surfaces Q , 0] as a possibility to provide an information 
about the basic processes of energy losses on microscopic 
level, important for the purposes of optimization in many 
technological applications. This is also connected with 
the search for best coating and lubricant materials. The 
principal physical picture, usually considered in relation 
with boundary friction on few molecular layers of lubri- 
cant liquid, is the sequence of " freezing- melting" pro- 
cesses on the lubricant, giving rise to discontinious (stick- 
and-slip) displacement of sliding surfaces 0] . Recently a 
new theoretical approach was proposed for microscopic 
sliding processes at extremely low velocities of motion 
and upon a monolayer of lubricant atoms 0, as can 
be the case for friction force microscopy (FFM). Based 
on the adiabatic formation of metastable states (simi- 
lar to dislocations in usual deformed crystals or defects 
in the Frenkel-Kontorova model and the ir follow- 

ing relaxation, this treatment shows how the (average) 
microscopic friction coefficient depends upon the mate- 
rial parameters of the contacting solids and lubricants 
and also how stick-and-slip jumps with atomic period- 
icity can develop in the microscopic friction force. It 
should be noted that a similar microscopic mechanism 
of dislocation-assisted sliding was recently proposed for 
contacting asperities in dry friction 0. 

Ultimately, with the discovery of almost spheroidal 
molecule of fullerene Ceo 01 (and/or cylindrical carbon 
nanotubes Q) a hope had arosen to use such closed 
molecular structures, as "free rotating" lubricants, for 
considerable reduction of the friction coefficient. For in- 
stance, a reduction of sliding friction coefficient was al- 
ready discussed, due to involvement of spinning motion 
of surfaces in contact [l0(. However, the FFM experi- 
ments with use of Ceo monolayers deposited over atom- 
ically smooth solid surfaces brought some contradictory 
results 0| ) To get their better understanding, a 

further theoretical insight is desirable on the elementary 



processes of boundary friction. 

The present communication is aimed on extention of 
the above mentioned adiabatic approach to the processes 
of boundary friction which include the internal degrees 
of freedom of the lubricant molecules. Within a simplest 
model, we make an attempt to show that, due to the dis- 
crete atomic structure of such a molecule, there are pos- 
sible qualitatively different modes of slow motion, either 
dissipative or non-dissipative, depending on the applied 
external load on the contact. For dissipative modes, there 
are energy losses resulting from stick-and-slip discontinu- 
ities, but these losses turn out much lower then for similar 
processes at sliding solid surfaces upon atomic lubricants. 
Besides, the mean value of friction force is found to be a 
non-monotonic function of the external load. 



THE MODEL 

Let us consider a two-dimensional model for the 
boundary friction system which includes two semi-infinite 
atomic arrays, the "solids" , with identical triangular lat- 
tice structure and a spatial separation d between their 
srufaces, and clusters of four atoms, the "molecular lu- 
bricants" , confined between the solids (Fig. ^) . The dis- 
tances between nearest neighbour atoms in both solids 
and in clusters are supposed invariable |l4| (that means, 
corresponding to "absolutely rigid" bonds), and the bond 
length for the molecule equals to the lattice parameter a 
for the solid surfaces. 

The model includes certain "soft links", which deter- 
mine the energy variation vs relative displacement of the 
solids. In the first turn, this is weak Lennard- Jones inter- 
action between an atom of solid and a lubricant molecule 
at distance r: 



V LJ {r) = e 
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characterized by the adhesion energy eo and equilibrium 
distance ro- Besides, a lubricant molecule can experience 
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FIG. 1: Schematic of the friction system: two solid arrays 
(grey circles) with identical lattice structures are separated by 
a four-atomic lubricant molecule (white circles) . The molecule 
center rests at the origin, and its rhombic configuration is 
determined by the angles 9 and <j>, while the symmetric dis- 
placements of solids are described by the parameter p (see the 
text). 



rhombic deformation which changes its energy as k sin 2 </>, 
where k and </> are the elastic constant and rhombic an- 
gle, respectively. We neglect the interaction between 
different lubricants, that is omit their possible collec- 
tive modes. Thus the molecular lubricants are supposed 
to contribute independently into the total friction force. 
Then we choose one particular molecule and set its center 
of mass as the origin of reference frame, so its configu- 
ration is determined by the "internal" rhombic angle, <f>, 
and the "external" orientation angle, 9 (say, with respect 
to the normal to interface). In such a frame, the two 
solids are supposed to be displaced symmetrically with 
respect to the lubricant molecule and their configuration 
is fully defined by the displacement parameter p chosen, 
e.g., as the smallest positive longitudinal coordinate of 
the atoms of upper solid. 

All this permits to write the full energy (per one lubri- 
cant) in the simple form: 



L-j(\r„ 
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k sin 



(2) 



where r n and R TO are respectively positions of atoms of 
solids and atoms in the lubricant molecule. Obviously, 
this extremely simplified model does not pretend to give 
a quantitative description and explanation of friction ex- 
periments in real systems with complex and typically in- 
commensurate solid and molecular structures. It serves 
mainly to illustrate some new qualitative possibilities for 
weakly dissipative processes, associated with the inter- 
nal (rolling or rotational) degrees of freedom of molecular 
(nanotubes including) lubricants, in contrast to the dis- 
sipation by only translational motions of the " point-like" 
(atomic) lubricants. 



THE ADIABATIC DYNAMICS 

The adiabatic treatment of the system, corresponding 
to the expression (J2J), follows the lines suggested in Ref. 
0. The equilibrium distance for Lennard-Jonnes inter- 
action Q is taken equal to the interatomic distance in 
the solids: ro = a. At given separation d between the 
surfaces, we calculate numerically the total energy pro- 
file E(p, 9, <j)) as a function of the displacement parameter 
p, and also of the angles 9 and <fi. Next this function is 
optimized with respect to the deformation angle (/) to re- 
sult in the profiles E(j>, 9), such as displayed in Figs. 12151 
(for different values of p, through the whole displacement 
period from to a, and at different separations d). 

The primary optimization in <fi refers to the stronger 
elastic deformation constant (we took k = 0.5eoj while 
the amplitude of relevant energy oscillations in 9 is 
~ 0.1eo)j an d hence to a faster relaxation in <\> than in 
fl.[l5j| Then the system behavior at very slow uniform 
variation of the parameter p with time (that is, the slow 
dynamics) is obtained from the analysis of the profile 
E(p, 9). Below we analyze how this profile changes with 
growing external load, which is here simulated by a grad- 
ual decrease of the separation distance d. 
i) At greatest separations (evidently corresponding to 
the lowest loads), it is seen from Fig. that the en- 
ergy profile has a single n/2 periodicity in 9. This means 
that, for any given p, there is a single equilibrium state 
for the molecule, characterized by its orientation 9 (p) 
(within to a C4 rotation) and energy E (9 (p)) , such that 
(dE/dd) e(p) = 0, (d 2 E/d9 2 ) g(p) > 0. In this case, the 
phase trajectory E(9) (shown by the sequence of bold 
dots linked by arrows, for growing displacement) is closed 
and continious. Thus, the system energy changes in a 
fully reversible way and, though some forces are exerted 
in the process, their mean value over the cycle and so 
the mean friction force are exactly zero. This reversible 
variation of the angle 9 around its median position (such 
that the long diagonal of the rhomb points vertically, in- 
set to Fig. [21, corresponds to a "rocking mode" of the 
molecular motion. 

it) When the two solid surfaces are getting closer (at in- 
creasing load), the first qualitative change that appears in 
the system behavior is the doubling of its energy minima. 
Now there are two splitted minima at the initial config- 
uration p — (Fig. |3J), which correspond to a two- fold 
degenerate equilibrium state of the deformed lubricant 
molecule: the long diagonal of the rhomb can deviate by 
a finite angle < ir/2 to both sides from the vertical (in- 
set to Fig. Since the molecule is considered a classic 
object, it initially occupies only one of the minima (the 
left one is chosen in Fig. |3J). 

However the degeneracy gets lifted for p > 0, so that 
one of the splitted minima turns stable (s) and another 
a metastable (m) equilibrium state. As it is seen from 
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FIG. 2: Series of profiles of potential energy vs orientation 
(at optimal deformation, see the text) and the related equi- 
librium states (bold circles) of a lubricant molecule at growing 
displacement of the solids (the numbers indicate p/a values) 
and at highest separation between them. The phase trajec- 
tory (arrows) shows that the molecule performs continuous 
"rocking" motions and after the full cycle (at p/a — 1) re- 
turns to its initial state with no energy dissipation. Inset: 
the rhombic configuration at p = 0. 
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FIG. 3: The same as in Fig. |2|but at closer separation between 
the solids (higher external load). The lubricant molecule after 
"rocking" returns to its initial state, but in this course it 
experiences a thermally activated jump (the dotted arrow over 
the potential barrier for p/a = 0.4) between the metastable 
(0.4) and stable (0.4') minima, and the corresponding energy 
difference is irreversibly lost. Inset: the configuration at p = 
corresponds to one of the splitted minima 



the consecutive curves in Fig. |2| the energy barrier h be- 
tween the m-state and the nearest s-state (here that to 
the left from m) decreases by many times with growing 
p, as a certain function h (p) . Since the adiabatic lifetime 
for the m-state is: r m = r a exp(/3/i) (where r a ~ 1CP 12 s 
is the atomic oscillation time and j3 is the inverse temper- 
ature) , it decreases in this course by many orders of mag- 
nitude. Eventually, this lifetime gets comparable to the 
characteristic time To of slow displacement by an atomic 
period, at a very sharply defined instant when p = po 
so that h(po) = ho — In (ro/r a ). Hence it is almost 
exactly at this instant that a thermally activated jump 
from m- to s-state is realized. After the jump, the en- 
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FIG. 4: The same as in Figs. [5]and[3]but at still higher load. 
The thermally activated jump changes its direction and the 
lubricant molecule rotates by the angle tt/2 after every period 
in p (" rolling mode" ) . The energy loss and hence the friction 
force is somewhat larger than that in Fig. |3] 



ergy difference between m- and s-states (to the moment 
of transition) is irreversibly lost, through the creation of 
quasiparticles (phonons, for insulating solids) which are 
finally thermalized in the bulk. For typical displacement 
rates in FFM ~ 10 2 pA/s, one has r ~ 10 -2 s, so that 
the barrier to the transition moment is still as high as 
ho ~ 23/3~ 1 . For the situation presented in Fig. 0] the 
transition corresponds to po = 0.4a (the potential profile 
shown by the solid line). It is also seen that, with fur- 
ther growing p, up to a, the system returns to its initial 
state; thus the phase trajectory is still closed ("rocking 
mode" again) though discontinious. This discontiniuity 
produces stick- and-slip features in the microscopic be- 
havior of the force vs displacement and it is the only 
source of the irreversible losses in our model system, in a 
full similarity to the model of sliding solids upon atomic 
lubricants. 

There are however some distinctions between the two sys- 
tems. Firstly, the stick-and-slip profile for molecular lu- 
bricants is more complicated than the simple triangular 
sawtooth for atomic lubricants. But especially important 
is the fact that the irreversible forces for molecular lubri- 
cants are smaller than elastic, reversible forces (the jump 
heights in Figs. 13151 are noticeably lower then the ampli- 
tudes of smooth oscillations), while the irreversible forces 
for atomic lubricants are orders of magnitude higher than 
the reversible ones. Since the reversible forces in both 
cases are characterized by the same energy scale £rj, this 
indicates a possibility to essentualy reduce the dissipa- 
tion by molecular (rotating) lubricants vs that by atomic 
lubricants. Note, at least, that if the molecule occurs ini- 
tially (at p = 0) in the right minimum, it simply spends 
a period, < p < a, in the s-state, but then passes to 
the m-state during the next period and eventually comes 
to the same closed discontinuous rocking. 
Hi) Now let us bring the surfaces yet a little bit closer, 
as by a very small change of the distance d between the 
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abrupt transitions from one regime to another at contin- 
ious variation of external load. In the sequence of regimes 
("rocking" and "rolling"), the obtained friction force vs 
load displays non-monotonic and hence non-linear behav- 
ior. This model (of course, with due improvements to be 
inserted) may provide a mechanism for explaining the 
data of real experiments with the fullerenes Ceo a s lubri- 
cants and give indications for an optimal regime of their 
practical applications. 



FIG. 5: The same as in Figs. EH at futher growing load. The 
molecule returns from "rolling" to "rocking" regime with an 
intermediate jump, but at much smaller energy loss than in 
Figs. OH and [I 

solids (Figs. |21 and 0J. Then the system evolution ac- 
quires yet a new quality: now at the transition moment 
the nearest s-state is to the right from the given m-state. 
This results in that after the transition, at p > po, the 
molecule does not return to its initial state, but will be 
eventually rotated by the angle 7r/2 (the dashed circles 
and arrow in the inset to Fig. 0J, and this rotation will 
be repeated each next period. Thus, in this case we have 
both unclosed (corresponding to the " rolling mode" ) and 
discontinious regime of lubricant motion. 
iv) With still increasing loads, the situation of a single 
minimum for E(0, 9) will be restored again (though dis- 
placed by the angle 7r/4), corresponding to a "square" 
molecule (inset to Fig. |3J. But now, unlike the case of 
low loads, this minimum gets splitted with increasing dis- 
placement p (here atpw 0.3). The molecule returns from 
" rolling" to " rocking" regime (" rock-n-roll dance" ) , but 
with considerably reduced energy dissipation: the energy 
loss at a jump for d = 2.1a (Fig. |3J) is about 2.5 times 
smaller than for d — 2.28a (Fig. i.e. it decreases with 
load increasing. We do not present here detailed results 
of numerical simulations for even higher loads, resuming 
only that they reveal a number of subsequent dissipative 
regimes, either with growing and falling friction forces. 

CONCLUSIONS 

The above simple analysis demonstrates that bound- 
ary friction with participation of molecular (spheroidal 
or cylindrical) lubricants can possess quite unusual prop- 
erties, such as existence of various regimes of molecu- 
lar motion, either non-dissipative and dissipative, with 
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